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Cold atoms in optical lattices allow to study in a detailed way a many body dynamics. Rapid 
time-dependent modifications of optical lattice potentials may result in significant excitations in 
atomic systems. The dynamics in such a case is frequently quite incompletely described by standard 
applications of tight-binding models (such as e.g. Bose-Hubbard model or its extensions) that 
typically neglect the effect of the dynamics on the transformation between the real space and the 
tight-binding basis. We develop a quasi-exact multi-band approach using time-dependent Wannier 
functions and apply it to model cases, directly related to experiments. 

PACS numbers: 67.85.Hj, 03.75.Kk, 03.75.Lm 
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Ultra-cold quantum gases in optical lattice potentials 
allow for precise studies of standard models known from 
other branches of physics (e.g. the condensed matter the- 
ory) as well as proposing novel situations with intrigu- 
ing properties. The latter utilize rich atomic internal 
structures, a versatility and an extreme controllability 
of atomic systems [lH3j- Many body physics often ad- 
dresses stationary properties such as phase diagrams - 
cold atoms enable also a controlled study of dynamics. 
That is especially interesting in the vicinity of quantum 
phase transitions Q where one of the intriguing problems 
is the adiabaticity and quantitative analysis of deviation 
from it for slow quenches (BHD]. 

Other interesting aspects of the dynamics concern ef- 
fects resulting from rapid changes of system parameters. 
The typical example of such a situation is a well known 
revival experiment [jj [To| where the sample prepared in 
a superfluid state is placed in the insulating environment 
by a fast increase of optical potential depth. Other pos- 
sible examples include fast quenches (Tl - 13 1 and periodic 
modulations of optical lattice potentials often used either 
for measuring the state of the system [H , 15 1 or even for 
modifying its effective parameters and thus its proper- 
ties |16H20|. Even faster modulations were suggested to 
efficiently populate excited bands (2l| . 

The aim of this letter is to show that for rapid mod- 
ifications of the optical lattice potential (e.g. its depth) 
a standard application of tight-binding models is incom- 
plete. We develop a quasi-exact multi-band theory which 
uses time-dependent Wannier functions and show its ap- 
plicability on few chosen model examples. 

For weak optical potentials, in the deep superfluid 
regime, the standard approach is to use a Gross-Pitaevski 
mean- field approach |22j,|23(. In deeper lattices, the de- 
pletion of the condensate becomes significant and an- 
other approach is necessary. The seminal work [24| uses 
a tight-binding approach mapping the real space system 
onto the lattice, resulting in a famous Hubbard model 



for fermions [25j or the so-called Bose-Hubbard model for 
bosons. One may consider 3D cubic lattices realized by 
three orthogonally polarized laser standing waves. Of- 
ten the reduced, two- (2D) and one-dimensional (ID) 
geometries are interesting flil . [l5| for which very deep 
lattices in remaining directions cut the atomic sample 
into 2D slices or ID tubes (with the confined degree(s) of 
freedom effectively described by the harmonic oscillator 
ground state). Explicitly, for the simplest quasi-lD situ- 
ation [H| Vi(f) = ssm 2 (kx) + ^-u! 2 (y 2 + z 2 ), Parameters 
s,lu are tunable in the experiment. 

Cold interacting Bose gas described by a second quan- 
tized Hamiltonian: 
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V 2 + V 1 (f)) 



+ iy dVdV* t (r)f t (r')l/(f,r')f(r)$(r'), (1) 

where V{r, r 1 ) is an isotropic short-range pseudopotential 
modelling s-wave interactions Q 
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with a being the scattering length. 

The field, *ff(r), is expanded in basis functions, 
W° L (f,s) = wf(x,s)H(y)H(z) [built as a product of 
Wannier functions in the direction of the lattice with 
harmonic oscillator functions in transverse direction: 
H{z) = (k/tt) 1 / 4 exp(-Kz 2 /2)[, n = mu/h: 



(3) 



i numbers the sites and a Bloch bands of the lattice. 

Performing the integrations in ([T]) using orthogonality 
of Wannier functions the extended Bose -Hubbard (EBH) 
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Hamiltonian is obtained 
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with J terms describing tunneling between sites while 
U terms 2-body collisions. Explicitly, Jf_^ = 

JW i «(r)(-£v 2 + F 1 (r))M7(f)d 3 f ! while Uffi* = 

Jd 3 rd 3 r"W l a (r^W^(f)V(r,f)W^(r)W l s (f'). The on- 
site energies Ef — Jq do not depend on site for trans- 
lationally invariant systems (we leave this dependence to 
signify that one may incorporate easily additional slowly 
varying inhomogeneous term). 

The tight binding representation is exact. The approx- 
imations emerge when we limit the number of bands as 
well as put restrictions on the Hamiltonian parameters. 
For sufficiently deep lattices (say s > 3) we may re- 
strict the tunneling to nearest neighbors only (see (27| 
for a shallow lattice case when next nearest neighbor 
tunnelings also play a role). Consistently, for interac- 
tions, we include terms such that i = j = k = I or 
i = j = k, and I being nearest neighbor of i up to 
a permutation (the so called density dependent tunnel- 
ings (2814301 are taken into account). From now on by 
EBH we shall denote the Hamiltonian @ with the fi- 
nite, low number of bands: a = 0, l,..,B. The corre- 
sponding Wannier functions are smooth and the action 
of the pseudo-potential is equivalent to a standard con- 
tact term V(r, f) = ^^5{r - f) = gS{f - ?). Re- 
stricting a to the lowest band only (a = 0) and taking 
solely i = j = k = I on-site interactions gives the stan- 
dard Bose- Hubbard (BH) model [24] . In the following we 
adopt the recoil energy En = k — 3? as an energy 
unit (A is a wave of the laser). We take fc _1 = A/2-7T as 
the unit of length. 

It is vital to note that Wannier functions depend on 
the lattice parameters, in particular s. While EBH gen- 
uinely describes the dynamics in such a lattice, we run 
into problems when, e.g. the lattice depth s varies in 
time. There are two options: either we keep the basis 
fixed in time determining it once at a given, say initial, 
s = so value or we make the basis time dependent so that 
Wannier functions change with time [i.e. we use s(t) in- 
stead of so m Eq. (jS}]. The former, while conceptually 
simpler, leads to difficulties: once s in the Hamiltonian 
(fT]) is different than so chosen for Wannier functions, the 
resulting Hamiltonian is no longer in the form of EBH as 
defined in (|4]) - in particular tunneling-type terms appear 
between different bands. Therefore most of the authors 
use the latter approach (at least implicitly) using the 
EBH (or BH just for the lowest band). Then changes of, 
e.g., lattice depth s(t) are just translated into changes 
of Hamiltonian parameters J"_^ and U"^ S evaluated for 



s = s(t) (see e.g. @, H, [HI, (H). 

Such an approach neglects the time-dependence of 
Wannier functions. The situation is similar to a ba- 
sic textbook unitary transformation case. Recall that 
if ip(t) = U(t)x{t) and the evolution of \ i s governed 
by the Hamiltonian H then the proper Hamiltonian for 
time evolution of ip is UHU^{t) + ihd/dt(U(t))W (t). Us- 
ing time dependent Wannier functions is equivalent to 
performing a similar transformation on our system. A 
straightforward calculation [26| yields the proper Hamil- 
tonian Hyv m the instantaneous Wannier basis in the 
form 

h w = n + w = n + ih^-(u(t))u f (t). (5) 

at 

Specifying further on that we are interested in changes 
of lattice depth in time we may write d/dtU(t) — 
d s {U)ds/dt and W = T ds/dt. Then we obtain [Hi 



t = -»e ir4-( s )K)S- 



T i-A s ) 



w"(x,s)—Wj(x,s)dx. (6) 



Transition integrals T"_j(s) obey relations: Vi,j,a 

a/3 _ ^_^ a +j3rpaP 



i—j 



0, TjZj = (-ir+PT^ = -Tp, T " p = for 
a — (3 odd. In particular, T"-term correction to a single 
band model (e.g. a standard BH) is zero. Thus the in- 
fluence of the W term is expected only when coupling to 
higher bands is appreciable. The T term contains both 
on- and off-site terms. Practically, most significant cou- 
pling occurs on-site between bands a — and a = 2. 
T^(s) decrease rapidly as |d| grows. In Fig. Q] most 
prominent of these parameters are shown as a function 
of s. In numerical examples below terms up to nearest 
neighbors are taken only. 

As mentioned above the W term is usually omitted in 
numerical simulations. Its relative importance depends 
on the product 7~4| . Clearly rapid changes of the lattice 
in time are necessary for the effects of W to be apprecia- 
ble. Thus for slow (e.g. 100ms) quasi-adiabatic quenches 
from shallow to deep optical lattices leading to Mott in- 
sulator phase formations 0, Q, [3, HH , the W term may 
be ignored safely. 

Situation becomes different for fast changes of s. We 
shall consider the influence of W term in model situa- 
tions. We shall discuss first a simple model of a linear 
quench. Then we briefly mention the famous revival ex- 
periment (Io| . Finally we analyze the recent proposition 
for efficient higher bands excitations [21j . 

A linear quench. It is realized assuming s(t) = sit/r + 
(1 — f/r)so, where r is the duration of the quench. Con- 
sider TV = 5 87 Rb atoms placed in a ID lattice of length 
L = 4 under periodic boundary conditions (PBC). The 
exact diagonalization gives the ground state at sq = 12 - 
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Figure 1: (color online) Relevance of different transition am- 
plitudes: (a) - nearest neighbor tunnelings J™ for different 
Bloch bands (b) - interaction integrals for g — l,n = 2n; 
the term Uuu° term present in the Bose-Hubbard Hamil- 
tonian is compared with interaction terms involving excited 
bands. Panels (c) and (d) show additional amplitudes [see 
Eq. ([§}] coming from the time-dependence of Wannier func- 
tions (i.e. the so called W term). 



the initial state. The quench is performed up to si = 40, 
for different values of the quench time, r. We find that 
(see Figure [2]) as soon as r < h/Eu the excitation en- 
ergy becomes significantly larger in the presence of the 
W term than without it. That reflects a significant differ- 
ences in the occupation of the second excited band. Thus 
a simple treatment of higher bands via the EBH model 
is insufficient to explain the dynamics; time- variation of 
Wannier functions has to be taken into account. 
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Figure 2: (color online) AE, the energy gain (i.e. the excess 
energy over the corresponding ground state) during a linear 
quench of a model fD system from sq = V2Er to s\ = 40Er,. 
For the adiabatic process AE = 0. Red (black, dashed) lines 
correspond to the simulation with (without) the W term in 
time-dependent Hamiltonian. Panels b) and c) show time 
variation of averages of band occupation operators {hi}, {hi}. 

In the limit r —> the quench becomes instantaneous, 
the evolution does not change the wavefunction. Yet the 
Wanner function basis changes from B(sq) to B(si) and 
so does the field operator, ([3]) representation in Wan- 
nier function basis. The basis change is realized via 



U{si)W (so) operation obtained directly from Eq.©. 

In the revival experiment [Io| a rapid quench is real- 
ized for bosons in the optical lattice by a rapid increase 
of the lattice depth. The authors were of course aware of 
the fact that too fast a quench would have populated ex- 
cited bands thus they chose the duration of the quench, 
t = 50/is, sufficiently long so that the higher bands ex- 
citations were negligible. As it turns out already for 
t = 20/is effects due to Wannier function dynamics are 
appreciable [26l |. 

Higher bands excitations. Recently Sowihski suggested 
[2~0 | to use periodic modulations of the lattice depth, say 
in the x-direction, 



:(*) 



j sin(o;£), 



(7) 



for excited orbital quantum state preparation. Large en- 
ergy gap between Bloch bands requires high frequency 
driving to couple different Bloch band states [2l|. This 
translates into significant values of ds/dt and contribu- 
tions from the W term in the dynamics may b e sig nificant 
- these terms were not taken into account in [21|. 

To see how the important the W par t is we have re- 
calculated the numerical simulation [21| using (Q| with 
and without the W term, Eq. ©. The studied system 
is a 2D lattice with the tight harmonic confinement in 
the third direction (26j. The system is assumed to be in 
a deep Mott insulator regime (s x = 32, s y = 20, k = 8). 
Due to a deep lattice potential, inter-site hopping can 
be neglected, and the whole system decouples into in- 
dependent 2D sites with the (assumed) integer filling 
v = 2. We have prepared the system in the ground 
state \ip{0)} with energy Eq. Following [2l[ we restrict 
the numerical simulation to first three bands (while this 
may not be sufficient for a simulation of a real situa- 
tion since W terms efficiently populate higher bands we 
consider the same model as in |2l| to isolate the influ- 
ence of Wannier functions' dynamics). We have per- 
formed numerical evolution of a system for a time cor- 
responding to 10ms with varying frequencies u x . As in 
[2l| we measure the maximal ground state depletion: 
8{u) x ) = l-sup te[0T ] \{i>uj :l {t)\'ip{0)}\ as a function of driv- 
ing frequency uj x . We find that the presence of W term 
changes significantly the depletion function in the fre- 
quency range considered (compare FigEJ) - Th e VV term 
leads to several additional excitations accompanied also 
by broadening and shifting the excitation peaks obtained 
without the W term. 

A key result of [2l| is a possibility of efficient popula- 
tion of higher Bloch bands. Population of excited state 
was found to oscillate with period of a few ms, making 
the whole process quite efficient. We have found that 
including the W term in the analysis makes populating 
the excited bands even faster. The oscillation period is 
decreased usually several times with similar excitation 
efficiency. Therefore, while confirming the possibility of 
direct resonant transfer of population to excited bands by 
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Figure 3: (color online) Excitation via modulation of the 
lattice height with and without W term. The depletion func- 
tion during the first 10ms, without (black, solid) and with the 
W term (red, dashed). The broadening of the peak around 
ui — 18.5 is a power broadening effect, see discussion in the 
text. 



lattice depth modulation, our analysis suggests that tak- 
ing the time variation of Wannier functions into account 
is crucial for controlling the process and for selective ex- 
citation of desired bands. 
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Figure 4: (color online) Floquet spectrum without W contri- 
bution (left) and with quasi-exact study, with W term (right 
panel). Broad avoided crossings for the latter are due to the 
strength of terms omitted within BH description as well as the 
influence of higher harmonics - see text for discussion. The 
region of avoided crossings (large curvatures) is highlighted 
with red (slightly lighter) color. 



The presence of the W term in the Hamiltonian may 
be also analyzed using the Floquet approach [16j . Ex- 
emplary spectra inspecting the broad structure around 
lu = 18.5 in the depletion function are shown in Fig. 2) 
The spectrum without W contribution shows a single iso- 
lated avoided crossing indicating a simple resonance (cor- 
responding to the isolated peak in Fig. [3]) in a contrast to 
broad structure of avoided crossings for the quasi-exact 
evolution. That structure correlates well with the broad- 
ened peaks observed in the depletion function. 

To understand this striking difference it is enough to 
consider the magnitude of different oscillatory terms in 
the Hamiltonian. For a deep lattice, when tunneling in 
negligible and a single site is considered only [2l| the 
driving comes from time dependence interaction term as 
well as from W term. The dominant interaction term 



C7 0000 may be expanded as [compare ([7])] 

QrrOOOO 

U°ooo {s) w ^oooo^O)^ SmSinu; t « 1.45+0.050 sinarf 



ds x 

for s° = 32 and s m — 4 (in recoil units). The dominant 
intrasite hopping term coming from the W term 



Sm sinwt s m ui cosuit 



dt V~ U " X/ ' ds, 

0.79 coswi - 0.026 sin 2ujt 



for lu = 18.5. Thus the driving term coming from 
W contribution is an order of magnitude stronger than 
driving induced by a modulation of interactions! The 
present of an "in phase" and "in quadrature" driving 
breaks the time-reversal invariance [34l| th at strongly af- 
fects the structure of avoided crossings |35ll36l |. Moreover 
W term brings strong second harmonic which addition- 
ally strongly modifies Floquet spectrum as well as the 
dynamics. For parameters of 2l| the W term dominates 
the dynamics. 

Let us mention also that the importance of the W term 
in the evolution appears also for sufficiently fast oscilla- 
tory movement of the lattice sites with a fixed lattice 
depth 



i6, m 



The selection rules for excitations are 
then modified with respect to the simple, translationally 
invariant, symmetric situation discussed here [26^ . Other 
possible applications of the presented formalism may be 
thought of due to the generality of the arguments leading 
to Eq. (|5|). In particular, with appropriate modifications 
it can be formulated for fermions as well. 

In conclusion we have shown that great care has to be 
taken when using discretized tight-binding approxima- 
tions for atomic systems in optical lattices in the presence 
of fast time variations of optical lattice parameters. A 
transformation to Wannier functions basis involves time 
dependent terms reflecting the time dependence of Wan- 
nier functions themselves. These terms modify the tight- 
binding Hamiltonian and are responsible for appreciable 
and directly experimentally measurable effects whenever 
the change of lattice parameters is sufficiently fast. 

It is worth mentioning that in a series of pa- 
pers, Heidelberg group [38l - l40j developed alternative 
approach for studying non-equilibrium dynamics using 
time-dependent variational basis which importantly takes 
interactions into account. The detailed comparison of 
both approaches is left for a subsequent study. 

We acknowledge discussions with Dominique De- 
lande and Krzysztof Sacha on the final form of 
this work. Projects of M.L. and J.Z. have been 
financed by Polish National Science Center un- 
der contract DEC-2011/01/N/ST2/02549 and DEC- 
2012 /04 / A/ST2 /00088 respectively. 
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SUPPLEMENTARY MATERIAL TO: FAST 
DYNAMICS FOR ATOMS IN OPTICAL 
LATTICES 

Mateusz Lacki and Jakub Zakrzewski 

Abstract: 

We provide a detailed derivation of the tight binding 
Hamiltoman which takes into account time-dependence 
of Wannier functions. We provide expressions valid for 
different dimensionalities of the optical lattice potentials. 



Multiband Bose-Hubbard Hamiltonian 

To make this material self-contained we repeat below 
some of the formulae from the letter extending it at the 
same time by pedagogical comments. 

Cold interacting Bose gas in an external potential V(r) 
may be described by a following second quantized Hamil- 
tonian: 



d 3 f^ t (f 



2m 



V 2 + V(f) *(r) 



+ - / d 3 fdV* t (r)* t (r')y(r,f')*(r)*(f') ) (9) 



H 



where V(f, f") is an interaction potential between bosons. 
For a very low temperature, the s-wave scattering dom- 
inates and the real potential may be represented as an 
isotropic short-range pseudopotential Q 



iitfra 



5{r 







d\r — F| 



(10) 



with a being the s-wave scattering length. 

We shall consider an external potential corresponding 
to a cubic optical lattice (generalizations to other lattice 
configurations are quite straightforward): 



V(r) 



^ s t sm 2 (ke t ■ r) 



(11) 



t=x,y,. 



where s t is the depth of the periodic standing wave po- 
tential in the t -direction, e t - versor in the corresponding 
direction, while k = 2tt/\ the wavevector for light with 
wavelength A. Consider x only and recall that Wannier 
functions, wf(x,s), localized at lattice sites, are linear 
combinations of Bloch functions 0™, eigenfunctions of 
the single particle problem with potential ssin 2 (kx) and 
quasimomentum p G [— k, k) (4l| : 



w?{x,s) 



-tk J ^»d P - (12) 

— k 



Now the field operator, *&(r), is expanded in terms of 
three-dimensional (3D) Wannier functions W"(r, s) (for 



the potential depths s = (s x ,s y ,s z )) that are products 
of the corresponding one-dimensional Wannier functions 



*(fO = 5>?WT(r, S ), 



(13) 



with 



W?(s,r) =w?*(s*,x)w^( Sy ,y)w«*(s z ,z). (14) 

In the expressions above the multiindex a = (a x , a y ,a z ) 
numbers different Bloch bands of the lattice while af 
is an annihilation operator for bosons at site i and in 
a band. Index i numbering the sites is, also really a 
multiindex i = (i Xl iy,iz)- 

Using the orthogonality of Wannier functions and per- 
forming integrations in © the extended Bose -Hubbard 
(EBH) Hamiltonian is obtained 



t—x,y,z . t . a 
a afi"i5 ijkl 

The notation i —± j, following (2l| in the sum indicates 
summation over all sites i and j shifted in t direction by 
distance i — j. The J terms describe tunneling between 
sites 

J*. = J W?(r) (-|^V 2 + V(r^j (r)d 3 r, (16) 
while U terms 2-body collisions 



U 



ijkl 



d Vd 3 r 1 W? ( r ) Wf ( ¥ ) V ( r , r 1 ) W% ( f) W^r 1 ). 

(17) 

The on-site energies Ef = J$ do not depend on site 
for translationally invariant systems (we leave this de- 
pendence to signify that one may incorporate easily ad- 
ditional slowly varying inhomogeneous term). 

Formally this representation of the Hamiltonian is ex- 
act. The approximations emerge when we limit the num- 
ber of bands as well as put restrictions on the Hamilto- 
nian parameters. In particular, restricting the expansion 
to the lowest band only, assuming that the tunneling am- 
plitudes are non-zero to nearest neigbors while interac- 
tions are entirely on-site we recover the standard Bose- 
Hubbard Hamiltonian 24|. Since we will be interested 



in processes populating excited bands we shall take them 
into account, limiting ourselves to several of them in each 
direction. 

As long as the number of bands included is limited 
to a finite number, any function built as a linear com- 
bination of Wannier functions will be smooth. Then a 
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pseudopotential, Eq. (|10[) . may be replaced by a Fermi 
delta-potential 



v 7 (r, r ) = o(r — r ). 



(18) 



From now on we shall assume a convenient units in 
which energy is expressed in recoil units (with recoil en- 
ergy E R = = and ^ = 1/fc is the unit of 
length. The potential, (JTHJ is then 



V(f, f) = g5(r — r), 



(19) 



with dimensionless coupling constant g = 8nak. The in- 
teraction integrals separate into a product corresponding 
to each coordinate taking the form 



vgr = 9 n 



ijkl 



t—x,y,z 



gU / dtw^(t)w^(t)wf(t)w^(t). (20) 



Frequently deep optical lattices are used to to separate 
atomic cloud into parts, a single retroreflected beam in 
one direction cuts a cloud into 2D slices, two such perpen- 
dicular standing waves produce a set of very weakly cou- 
pled tubes. Separating different energy scales one may in 
those situations assume, that in the tightly confined di- 
rection (directions for tubes) the system remains in the 
lowest band and neglect the tunnelings between slices 
(tubes). This has to be done with care, virtual effect 
of high lying excited bands may be not negligible [28- 
|30| . In this limit one may often approximate the Wan- 
nier function in that lowest band, say in z direction by 
a ground state of an appropriate harmonic oscillator, i.e 
Wi" ~ H(z) = (k/tt) 1 / a cxp(— kz 2 /2), where dimension- 
less k = hui z /2En = \J Sz/Er where s z corresponds to 
the lattice depth in that direction. Thus the basis func- 
tions for 2D system with the effective potential of the 
form 



V-z{f) = s x sin 2 (fcr) + s y sin 2 (fcy) + y w z^ 2 

are 

W« (r, s) = <* (x, 8 x )w? (y, s y )H{z). 
For ID tube, the corresponding formulae read 



and 



Vi(r) = ssm 2 (kx) + y^ 2 (2/ 2 + z 2 ), 



W z a (f,s) = wf{x,s)H(y)H(z) 



(21) 



(22) 



(23) 



(24) 



Observe that in the reduced geometries tunneling takes 
place along the periodic lattice potential only. This 
simplifies the notation, in particular in ID geometry, 
the index a numbering the bands is simply an inte- 
ger, a = 0,1,.., B, the tunneling is unidimensional and 
Hamiltonian (4) of the paper is obtained as an extended 
Bose Hubbard system . 



Contribution of the time dependence of Wannier 
functions to the tight-binding Hamiltonian 

Suppose that we are interested in the dynamical prob- 
lem with time dependent lattice depth s = s(t),s(0) = 
so- The initial wavefunction is expressed in Wannier basis 
{W"(-, so)}. Similarly the EBH Hamiltonian is obtained 
for that value. As discussed in the letter if s changes 
in time so does both the tunneling and interactions pa- 
rameters of the EBH Hamiltonian as well as the Wannier 
functions themselves (an alternative approach of keep- 
ing the basis fixed in time would require a different tight 
basis Hamiltonian). Therefore, when s = s(i) is time de- 
pendent so is the isometric basis transformation lA(s(t)) 
from the position representation to the lattice (Wannier) 
representation. Of course H = U(s(t))Hx^ (s(t)). De- 
fine the transformation via ip(t) = U(s(t))ipx(t) (where 
ipx{t) is the wavefunction in the position representation 
while ipit) corresponds to the lattice). Then a standard 
textbook derivation gives the proper new Hamiltonian 
T-Lw in the form: 



n 



H + ih(J^U{s(t)))tf(s(t))=U + W, (25) 



and the TDSE for ip(t) : ihd t ip = U w i>- 

For s(t) changing in time slowly enough, the second 
term may be neglected (this corresponds physically to 
the assumption that the system has time to adapt to a 
given change of basis). For quick variations of s(t) this 
term leads to appreciable effects, we shall now evaluate 
its form. 

The basis for the Hilbert space for a gas of N bosons 
consists of symmetrized (tensor) products of N single 
particle basis functions W". In restricted lattice geome- 
tries, we study of multiparticle states which in transverse 
direction(s) contain a harmonic oscillator ground states 
- we take W functions as in (|22|) or (|24|) . These wave- 
functions are effectively ID or 2D with finite, effective 
"width", set by the curvature of transverse confinement. 
We denote such a basis for N particle problem by W N (s), 
the basis depends on the lattice height by the value of 
s parameter through the set of single particle Wannier 
functions W(s) = {W?(; s)}i, Q - 

The Hilbert space for the lattice system (in which 
the Bose-Hubbard hamiltonian is usually expressed) has 
time-independent basis (the Fock basis) T . 

Let us define a shortened notation. The lattice Fock 
state with occupation n L of mode i will be denoted as \ ftc) 
The corresponding state in the position representation for 
the lattice with height s will be abbreviated: \nx , s). We 
always assume that m + . . . nc = N. 

Action of the map lA(s) from the continous space with 
base W N (s), to the Fock space with base T is rather 
trivial: in the chosen orthonormal bases it is the identity 
matrix — it maps state \rix,s(t)} to a state \nc,s(t)}. 
The map is thus always an isometry (note that for ID 
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and 2D lattice it is only a partial isometry from a full 
continuous space). 

Single particle states defined by Wannier functions and 
in the discrete lattice are enumerated by two indices a 
and i, from this point up to the end of the derivation, we 
introduce the multiindex i = (a, i) to simplify notation. 

Using this notation we can express the map IA as: 



U(s) = ^2 \ r nc){mx,s\ 



Now we expand the derivative of U isometry. We 
use the fact that basis of the Fock space T is time- 
independent: 



—i 

dt 



U(s(t))\ \i>) x =Y,\mc 



dt 



[m x ,s(t)\ |V> 



correspond to changing the mode of only one particle 
from configuration ft — the mode t to x. Therefore: 



rrii — rii i =/= i, k 

m t = n L — 1 



(29) 



As T L L = due to norm preservation, only t ^ ^ terms 
contribute. Change of occupation is compatible with ac- 
tion of al c a L operator. We will show that also the numer- 
ical factor agrees. A mode to be differentiated (mode t) 
may be chosen in (|27[) in n L ways, and: 

n L yJn L (n*c + 1) 



^/m!n 2 ! . ..n Ll \ ^Jn x \ . . . (n L - l)\ . . . (nZ + 1)1 . ..n Ll l 

(30) 

Thus from (l26l). (l27l) and (l28l): 



Thus the W term is just: 



W = ih^2\rft £ )(ft c \ (j f (m x ,s(t)\j \ft x ,s(t)) = 

n.m 

= -ihJ2\™c )(nc\(m x ,s(t)\ (^\ft x ,s(t)))(26) 



The relation between ft and rh that may give nonzero 
contribution to the above sum remains to be worked out 
as well as exact values of the coefficient. To do so, we 
expand the time derivative, by inserting exact action of 
the symmetrization operator: 



E 



dt 



^N\m\n 2 \ . . . n t0 !— \H X , «(*)> = 
at 

[Wi(x t( i ) )W / "i(x w(2 )) . . . Wl(Xn( ni )) 



■rrES(N) 

^7r(ni + l) 

)...w 2 ( )...w l 

\^tt(N) J J 



k=l TreS(N) 

■W 2 {x^ ni+r) ) . . . W L (x^ (k) ) . . . W Lo (x„ w ) (27) 

In the above line each of Wannier functions W L depends 
on t through s(t). The formula is well-stated, because 
only finite number of modes has nonzero occupation: for 
1 > io we have nj = 1 and no factors W L . 

Now we use the partition of unity \W X >< W x \, 
appying it on W L we get: 



-W L {x) = Y,T?W„{x) 



(28) 



where T* = J W lt (x)W L (x)dx. Therefore, by combining 
together (gSJ), (gT]) and igSJ) one obtains that the only 
(ft, fh) giving nonzero contribution in (|28[) are those that 



\ft x ,s(t)} = T?^/n L (n„ + l)\m x ,s(t)) (31) 



dt 



Above m is assumed to satisfy relations (|29|) . All in all, 
we obtain 



U(t)U\t) = -J2 T ^c)(ncWno(n tc + 1) = 

1 .x 



We now go back to the original labeling by Bloch band 
number: 1 = (a, i), k = (/3,j). To obtain the form of the 
VV term from the 

Now T* = if", = -T afl .. We obtain the form of W 
term used in the main article: 



w 



3^5 = / W?(x)W?(x)d*x. 



(33) 



The term T* has to be worked out for the basis 
functions for the lattice in the appropriate dimension. 
In the ID lattice, Wannier functions are of form (f24|) . 



then: T* = J dydzH(y) 2 H(z) 2 fdxw?(x)w? (x). 
Due to normalization: T* = J dxWj(x)wf(x). 
For 2D lattice, from ([221), we S et: T ? = 
JdzH(z) 2 Jdxdyw^(x)w^(y)f t (w£ (x)w% (y)) . 
Nonzero values may be obtained only if i x = j x Aa x = /3 X 



or 1. 



jy A Oty 



fly. Thus the W term per- 



form hopping of a particle in only one direction 
(including the Bloch band change). The corre- 
sponding amplitude for hopping in y direction are: 
J dzH(z) 2 J dxw™J (y) 2 f dyWj (y)w" y (y), which after 
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normalization becomes just: J dywj (y)w°' y (y). Simi- 
larly for the x direction we obtain: J dxiVj* (x)™"" (x) . 
Analogously for the 3D case fusing IT4"|) the amplitude for 
hopping in direction t £ {x, y, z} is: J dtwj* {t)ii\* (t). 

Analysis of the revival experiment 



20 50 [us] 1 2 3 4 5 




Time [ms] 



Figure 5: (color online) Time dependence of the order pa- 
rameter in the lowest band, <j> = (i!{t)\a°\%b(i)) for 50(is initial 
linear quench as used in the experiment flOl ] (top panel) and 
for a faster quench of 5/j.s (bottom panel). Red full line - ex- 
act evolution using Hw Hamiltonian, Eq. (|25|l . black (dashed) 
line the result obtained neglecting W term. The difference 
between the two is negligible for the top panel, showing that 
the higher bands play a minor role in the experiment. For the 
faster quench (bottom panel) the W term describing time- 
dependence of Wannier functions is crucial for the description 
of dynamics. Note that the time scale during the quench is 
expanded in comparison to the subsequent evolution to make 
the plot more readable. 

In the revival experiment flo| the atomic system is ini- 



tially prepared in a relatively low lattice, on the super- 
fluid side. Then a rapid quench is realized by a sudden 
increase of the lattice depth. The subsequent evolution 
of the system in a deep lattice is monitored by measur- 
ing the time dependence of the contrast of interference 
fringes in the momentum distribution. The initial coher- 
ent state like occupation of sites evolves differently in the 
deep insulating-type lattice showing the decay and par- 
tial revivals of coherence. In the simulation the evolution 
of coherence may be monitored by the time dependence 
of the order parameter <f> = (a ) where the superscript in- 
dicates that the lowest Bloch band is taken into account 
only (we drop the site index, as we shall consider a single 
site only - in a very deep lattice the tunneling may be 
neglected and a single site evolution is considered - see 

The contribution of W term (see the letter) depends on 
the speed of the quench, it becomes important when ex- 
citations of higher Bloch bands become appreciable. The 
authors (lo| wanted to avoid population of excited bands 
(that could be controlled experimentally) thus they ex- 
perimentally chose the duration of the quench, r = 50/is, 
sufficiently long so that the higher bands excitations were 
negligible. Consequently a contribution of Wigner func- 
tions dynamics (the W term) is for experimental param- 
eters quite small (see the upper panel in Fig. [5] As it 
turns out already for r = 20/is effects due to Wannier 
function dynamics are appreciable while a still faster 
quench leads to strong Wannier functions dynamics ef- 
fects (compare Fig. [5p). Let us point out that to make 
calculations less computer demanding we used a 2D lat- 
tice, the initial state, ip(0) was prepared as a coherent 
state with (no) = 2 at sq = 8 and a linear quench up to 
si =40 was realized. That roughly corresponds to the 
revival plot in Fig. 2 of [Icj . 

It seems, therefore, that making the quench time in 
the experiment [lo| shorter by an order of magnitude 
would allow for a direct experimental verification of the 
approach presented here. 



